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Abstract 

We present an enumeration of all possible amplitudes that contribute to an n-jet process 
in QCD. We estimate the number of amplitudes for large number of jets and determine the 
actual number of amplitudes to be calculated, which is smaller due to relabelling among 
(massless) quark flavours. 

1 Introduction 

With the advent of high-energy hadron colliders such as the Tevatron and the LHC, there arises a 
need for accurate QCD calculations of amplitudes of increasing complexity, either as problems in 
their own right, or as possible backgrounds to other physics. The complexity becomes apparent 
either in the number of loops to be considered, or in the number of external legs in the diagrams: 
the present paper aims to deal with the latter of these issues. In recent years there has been 
considerable progress in the computation of multi-leg QCD amplitudes [§, ^ Essentially 
based on the earlier work of [||, |6|, these new algorithms employ the recusrive structure of the 
Schwinger-Dyson equations to express the full amplitude in terms of smaller subamplitudes 
in the essentially most compact way, thereby reducing the computational complexity of these 
amplitudes from roughly k\ to about 3 fc , where k is the number of external legs. This enormous 
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improvement has led to the becoming feasible of amplitudes with as many as 8 or 9 outgoing 
partons. 

Another issue then arises, that of summing the contributions of all possible QCD processes 
to a given multi-jet final state. As we shall show, the number of QCD amplitudes contributing to 
the probability of an observed event increases very rapidly with the number of jets, so that the 
following questions become relevant. How many amplitudes contribute precisely? What is the 
asymptotic form of this number for large multiplicity? To what extent is the folk-lore that the 
purely gluonic amplitude dominates the cross section still valid? Given that many amplitudes 
can be related to each other by a simple relabelling of the quark flavours, how many distinct 
amplitudes have to be computed? These are the questions addressed in this paper. 



2 Physical complexity: contributing amplitudes 

Under the assumption that the various (anti)quark types and gluons cannot be distinguished 
experimentally, and that all these parton types are (essentially) massless, the only information 
experimentally available about any given event is the configuration of the observed momenta. 
We shall denote such an event by its momenta as follows: 

P1+P2 -> qi + Q2 + <?3 H 1- q n , 

where n outgoing partons/jets are observed. To obtain the total probability density for this 
event in phase space, one has to consider all possible 2 — > n QCD amplitudes[|,viewed here as 
functions of 2 + n momentum arguments, and assign the observed momenta to these arguments 
in all possible ways without double counting. Note that, due to the composite nature of the 
incoming hadrons, also the initial state may require more than one assignment: for instance, 
a quark-gluon initial state q(pi)g(p2) is to be counted as distinct from q{p2)g{p\), whereas of 
course the purely gluonic initial state g{p\)g{p2) is counted only once. Similarly, if the final 
states contains m quarks of a certain type, a corresponding factor 1/m! has to be applied. In 
what follows we shall denote the number of (essentially) massless quarks in the final state by /, 
so that / = 4 at relatively low momenta where the b quark might be identifiable, / = 5 typically 
for QCD studies at the LHC, and / would be 6 at some future multi- hundred TeV collider. The 
number of flavours contributing appreciably in the initial state is denoted by j, so that j = 3 
would be appropriate if the charm quark structure function can be neglected, and j = 4 if it is 
included. We shall, however, keep to general j and / as much as possible. 

2.1 Arrangement of the initial states and final states 

The various possibilities for the initial states are: 

x Of course, squared and spin/colour-summed and -averaged in the usual way. 
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• gg Obviously, there is only one possibility for the initial state. The final state can be 
anything so there are n partons in the final state. 

• qiqi If we have j flavours then there are 2j possibilities for the initial state and n partons 
to be distributed in the final state. 

• 9Qi ■> 9Qi There are 2j initial states as in the previous case. For the final state, we know 
that there must at least one qi or q\ so there are (n — 1) + qi or (n — 1) + qi partons. 

• QiQi i QiQi For the scattering of identical quarks (anti-quarks), there are j initial states and 
since the same partons must appear in the final state we can have n — 2 partons plus the 
initial quarks or anti-quarks. 

• QiQk i QiQk 5 i 7^ k For the scattering of different quarks (anti-quarks) we have j(j — 1) 
possibilities for the initial state, and again n — 2 partons plus the initial quarks (anti- 
quarks), for the final state. 

• QiQk ,i^k For this final case we have 2j(j — 1) initial states and n — 2 partons plus the 
quark and the anti-quark in the final state. 

All of the above can be summarized in the following table. 



Initial States 


# possibilities 


Final States 


gg 


1 


n 


qm 


2j 


n 


gqi 


2j 


(n - 1) + qi 


gqi 


2j 


(n - 1) + q~i 


qiqi 


3 


(n-2)+qi + qi 


qiqj,i^ j 


3(3 ~ 1) 


(n-2)+qi + qj 


Mi 


3 


(n-2)+qi + ^ 


Mj^ + 3 


3(3 ~ 1) 


(n-2)+q~i + qj 


qiq~j,i^ 3 


2j(j - 1) 


(n-2)+qi + qj 



where i,k = 1, . . . ,j. From the second column we can read off the total number of initial-state 
momentum configurations: 

1 + 3(2j) + j + j(j - 1) + j + j(j - 1) + 2j(j - 1) = (1 + 2jf (1) 

From this table we can arrange four different groups of initial states, which differ in the flavour 
structure of their final states. They are shown in the following table: 
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Group 


Initial state 


# of final states 


A 


99, Wi 


A(n) 


B 


gqi, 9^ 


B(n) 


C 


qm, Mi 


C{n) 


D 


mk, Mk, q-iq~k,i + k 


D{n) 



2.2 Counting of the final states 

Group A The distinct possibilities for flavourless final states are: 

n = n *(g)+ni* (qiqi) + n 2 * (92^2) H h n f * (g/g/) (2) 

where no is the number of gluons g and ni,rt2, . . . ,nj are the numbers of qf and q~f quarks 
with different flavour /. The number of different processes A(n) is the number of the various 
dinstinct ways to distribute n different final momenta among n partons: 

(n)! 



A(n) 



E 



no,ni,...,nf>0 



(n )\( ni )P(n 2 )l2...(n f ) 



[2 



G(n + 2ni + 2n 2 H h 2n/ = n) 



(3) 



where 0(a = b) = 5 a ^. We can evaluate this number by forming the generating function: 



= E ^(*) 



E 



x 2n f 



k>0 



no,ni,...,rif>0 



no } - (ni)! 2 (n 2 )! 2 (n/)! 



12 



K n>0 IL ' J \n>0 

where Iq(x) is the modified Bessel function of the first kind and zeroth order 




(4) 



Group B All the possible final states for this case have a single net flavour, and can be 
written as follows 



n = n *(g) +n 1 * (qiqi) +n 2 * (g 2 <?2) H h n f * (g/g/) + g» 

The number of different processes B(n) is: 

(n + 2ni + 2n 2 H + 2n/ + l)! 



(5) 



Bin) 



E 



no,ni,...,riy>0 

This gives the generating function 



B(x) 



n !(ni + l)!(ni)!(n 2 )! 2 ---(n/)! 2 



-e(n + 2ni + 2n 2 H h2n / + I=n) (6) 





/-1 



E 



„2n+l 



(n)!(n + l)! 



e x •I (2x) / - 1 -/ (2a 



(7) 



v n>0 / \n>0 v y / \n>0 

where the prime denotes the derivative of the Bessel function with respect to the argument 2x. 
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Group C For this case the final state is 

n = 2n * (g) + m * (qiqi) +n 2 * (<? 2 <?2) H h n/ * (<?/<?/) + 2 * (%) (8) 

and the number of possibilities is 

The generating function is 



C M = E — C W = e"/o(2x)^ 1 • {2/^(2x) - J (2x) (10) 

n>0 n - 



Group D The derivation goes through as in the previous cases and the result is 

V{x) = £ ^D(n) = e-/ (2x)/- 2 • (l>(2x)) 2 (11) 

n>0 n - 

The total number of possibilities for the final state can now be determined: 

G(n) = (1 + 2j)A(n) + AjB{n) + 2jC(n) + 4j(j - l)D{n) (12) 
with the generating function 

= E -T G (") = &X i (l + 2j)/o(2^) / + 4i/ (2x/- 1 ^(2x) 

n>0 n - 

+ 2j/ (2x)/- 1 (2^(2 a ;)-/o(2x)) 

+ Aj{j - 1) 7 (2x/- 2 (^(2:r)) 2 } (13) 

We can put this in a more compact form: 

G(x) = e x I (2x)^ (l + I (2xy (14) 

To get the number of processes we expand the generating function Q{x) and pick out the relevant 
coefficients. For example, for / = 3, 4, 5 flavours we have: 
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Total number of amplitudes 




/ = 3 


/ = 4 


n 


J = 2 


J = 3 


i = 2 


J = 3 


i = 4 


2 


71 


127 


81 


141 


217 


3 


299 


511 


377 


625 


921 


4 


1,763 


3,301 


2,645 


4,867 


7, 761 


5 


8,955 


16,297 


15,325 


27,087 


41,889 


6 


54, 353 


103,279 


113,733 


213,879 


345,465 


7 


304, 701 


570, 367 


745,421 


1,364,811 


2,162,617 


8 


1,879,723 


3,595,177 


5, 704, 061 


10,836,831 


17, 605,249 



Total number of amplitudes 






/ = 


= 5 




n 


J = 2 


j = 3 


j = 4 


j = 5 


2 


91 


155 


235 


331 


3 


455 


739 


1,071 


1,451 


4 


3,647 


6,601 


10,419 


15,101 


5 


23,255 


40, 157 


61,059 


85,961 


6 


200, 473 


372,719 


598,005 


876,331 


7 


1,470,061 


2,636,375 


4,118,865 


5,917,531 


8 


13,229,719 


24, 937, 645 


40,333,059 


59,415,961 



2.3 Gluonic contributions 

An interesting question that arises is the issue of contribution of gluonic processes, compared to 
the total number of processes, since often the purely gluonic process is assumed to be typical or 
'dominant'. In particular we would like to examine to what degree purely gluonic amplitudes 
dominate over other kinds of processes, since gluons have a different color charge than quarks 0. 
To this end we assign to each external gluon an additional factor k, resulting in a modification 
of the generating function ([14T): 

G k (x) = e kx h{2x)f-i {k + ^j 2 I (2xy (15) 

The corresponding generating function for purely gluonic processes would be: 

g° k (x) = k 2 e k * (16) 

We may compare the coefficients of the expansion of the two generating functions. This can 
be seen in the graph that follows, where we have plotted the ratio Q^, n jQ k ^ n of the coefficients, 

2 Note that we do not address the question of the singularity structure of gluonic versus other amplitudes. 
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for number of jets ranging from n = 2 up to n = 8, against k, and for the particular case of 
/ = 3, j = 3. Notice that the ratio approaches one as the factor k grows larger, but decreases 
with n. 




10 20 30 40 50 60 70 80 



In order to estimate how large k has to make the gluonic amplitude the dominant one, we look 
for values of k that give Q\ n IQ kn = 1/2. These can be seen in the table that follows, for the 
case / = 3, j = 3 again and for various numbers of jets. 



n 


2 


3 


4 


5 


6 


7 


8 


9 


k 


5.72 


7.13 


8.47 


9.78 


11.04 


12.27 


13.48 


14.67 



Another extension would be to account for the fact that the gluonic structure function is typically 
larger than that for a quark. Then the generating function becomes: 

G s (x) = e x I (2x)f-i (s + ^ J (2x)' (17) 

where S denotes the gluon structure function enhancement factor. We see that this is a mono- 
tonically increasing, quadratic function of S. For large S the generating function can be approx- 
imated by 

G s (x) = S 2 e x h{2xY (18) 

The coefficients of Gs{ x ) m ay be written Q n = C n S 2 where the C n depend on n. We can estimate 
the 'strength' of the 'S-extended' gluonic amplitudes compared to purely gluonic processes, 
Gn = S 2 , by computing the ratio r = lim s ^oo Gn/Gn- These ratios are tabulated below: 
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The ratio (%/Q n 


n 


/ = 3 


/ = 4 


/ = 5 


2 


0.1428 


0.1111 


0.0909 


3 


0.0526 


0.40 


0.0322 


4 


0.0078 


0.0046 


0.0030 


5 


0.0019 


0.00108 


0.00068 


6 


0.0003 


0.00012 


0.000066 


7 


0.000061 


0.000023 


0.000011 


8 


0.96810~ 5 


0.28910~ 5 


0.11410~ 5 



We conclude that for sizeable n the purely gluonic amplitudegives only a very small contribution. 
2.4 Asymptotic results 

It may be interesting to estimate the number of amplitudes for large number of jets. To this 
end, we would like to obtain the asymptotic form of the generating function for large n. The 
asymptotic expansion for Iq(2z) is 



a 2z 



I (2z) 



(2n)! 2 



z — ► oo. 



64™ n! 3 ' 

This expansion holds for Re(z) > 0, but we also have z) = Iq(z). For the function 

f(x) = I (2xf 

the asymptotic expansion is 

T ni T n 2 " " " T np @(ni H h n p 



n>0 X 



n 



(19) 



(20) 



(21) 



ni,...,ra p 



where N = (An) 2 . So the derivatives in the generating function Q(x), read: 
(l + /(*) = /(*) + 2/'(x) + /" (s) = Ne 2 ^x^ § 



where 



/3 n = (1 + 2p) 2 a„ - 2(1 + 2p)(n + I - l)a n _i + (n + | - l)(n + | - 2)a n „ 2 , 



and for the generating function we get 



g( x ) 



(4vr)// 2 



^(1+2/) 



C //2 



n>0 



2^ x n 



(22) 



(23) 



(24) 
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where 

In = @(ni+n 2 = n) a ni (f - j) /?„ 2 (j) (25) 

ni,n2>0 

The first few 7's for various numbers of initial and final flavours, are shown in the next table. 



Number of flavours 


7o 


7i 


72 


73 


/ = 3 


J = 2 


25 


-85/16 


249/512 


1873/8192 


i = 3 


49 


-189/16 


177/512 


1337/8192 


/ = 4 


J = 2 


25 


-15/4 


19/32 


101/256 


i = 3 


49 


-35/4 


15/32 


109/256 


j = 4 


81 


-63/4 


3/32 


93/256 


/ = 5 


3 = 2 


25 


-35/16 


409/512 


4871/8192 


i = 3 


49 


-91/16 


401/512 


6143/8192 


J = 4 


81 


-171/16 


273/512 


6831/8192 


j = 5 


121 


-275/16 


25/512 


6935/8192 



The most important term in the series is of course the first 

70 = a (/ - j) 0o(j) = a (f - j) (1 + 2jf a (j) = (1 + 2jf 
Keeping only this term in the generating function we have 



g{x) 



1 



(4tt)// 2 



x(l+2/) 



1 



rf/2 



(l + 2j) 2 , x^oo 



(26) 



(27) 



Let us, now, assume that we want to include the first K terms in the asymptotic expansion 
of Q(x), that is, we set jj to zero for j > K. The Borel transform 



t_XJ 

F{x) = j dyy K+f/2 e- y g(xy) 



has the expansion 



F(x) = T ( n + K + f/2 + l)g ri 

n>0 



X 



(28) 
(29) 



where G(x) = J2 n >o^nX n ; our approach consists in finding the coefficients of T{x) by studying 
its singularity structure. The integral for T(x) can be written as 



F{x) 



^Jdy exp[{-y + xy(2f + 
lk(K-k)\ 



(4tt)// 2 



lk- 



k=0 



,K-k 



•k+f/2 



K 

E 



(4vr)// 2 ^ x fc +// 2 (l - x(2/ + \))K-k+i ■ 



(30) 
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This expression is has a pole at xq = 1/(2/ + 1). Note that, due to our use of the factor y K+ $ / 2 , 
the integral (|2~8| ) is indeed dominated by large values of xy when x approaches x$, thus justifying 
the use of the asymptotic expression (p4|). Furthermore, there is of course a similar singularity 
which appears when we use negative x values: however, since that is located at —1/(2/ — 1) 
and hence further away from the origin than xq, this pole will give exponentially suppressed 
contributions which will not show up in our result for Q n . The k^ term in the series for J-(x) 
is seen to contain poles at x = xo of order up to and including K — k + 1: 

x -*-//2(l _ ±)-K+*-l = 
XQ 

The dominant behaviour of the coefficient of x n in the series expansion of this term is, therefore, 

1 ^ (k + f/2 + r)\ (n + K -fc-r)! _ 
^ r \(k + f/2- 1)! n\{K-k-r)\ ~ 

1 (n + K + f/2)\ 

x n+k+f/2( K _ k y( n + k + f/ 2 )\- [6 ' 

Inserting this in the expression for J-{x) and dividing by the factor T{n + K + f/2 + 1) to get 
the coefficient Q n , we see that K drops out from the expression, so that we may take it as large 
as we please. The resultant form for Q n is, therefore 

r r asy (2/ + l)"+// 2 ^ 7fc(2/ + l) fc 

In order to estimate how accurate this asymptotic expansion is, we have calculated the ratio 
between the exact and the "asymptotic" number of processes. The results are shown in the next 
table, where we have recorded the way these numbers improve as we add more terms in the 
asymptotic expansion of the generating function. Thus uq is such that G n /G^ sy is between 0.95 
and 1.05 for all n > no, when we include only the first term in the expansion, i.e. the term that 
contains 70. Similarily, n\ is the number of jets when we include 70 and 71, 712 when we include 
70, 71 and 72, etc. 



n 
x 
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Number of flavours 


n 


ni 


n 2 


ri3 


f Q 

J — 3 


J = 2 


26 


6 


6 


6 


J = 3 


31 


5 


4 


5 




J = 2 


21 


8 


8 


6 


/ = 4 


J = 3 


28 


8 


6 


6 




J = 4 


31 


7 


7 


7 




J = 2 


11 


12 


10 


8 


/ = 5 


J = 3 


21 


12 


10 


10 


J = 4 


25 


12 


10 


10 




j = 5 


29 


12 


11 


10 



Note that on some occasions, like for example / = 5, j = 2, the number increases as we add 
more terms in the expansion. But this is due to a small increase of the ratio which is improved 
immediately when we add the next term. 



3 Computational complexity: distinct amplitudes 

In the above we have shown how all amplitudes contributing to a certain cross section can be 
enumerated. This would also, then, be the computational complexity in an approach where each 
amplitude is calculated from scratch. However, there is of course a simplification owing to the 
fact that amplitudes that differ only by a relabelling of the (massless!) quark flavours are equal 
apart from a trivial difference in the structure function. It therefore behooves us to take this 
simplification into account. Now it must be kept in mind that, when a quark flavour occurs 
in the initial state, we should not relabel it since that is also taken care of by the factors 2j, 
j(j — 1) etcetera in table 1. Only those quark flavours that do not occur in the initial state may 
be relabelled. Let us perform the relabelling in such a way that the relabelled flavours occur in 
order of increasing multiplicity. As an example, in the process gg — ► X this means that instead 
of 

n \ x no+2(ni^ \-rif) 

^■^5SFW' (34> 



we have to determine, rather, 



n ] x no+2(niA hrif) 

Af{x) = £ E noKm!) 2 ---^ ; (35) 



no>0 0<ni<ri2<- ■ -<n f 

likewise, for the process gq — > X we have to compute 

n!x n +l+2(n 1 +...+n / ) 
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It is important to note that those quark flavours that can be relabelled occur symmetrically in 
these sums. 

The straightforward implementation of such inequalities appears to lead to horrendous com- 
plications. An exception is the following generating function: 

Z f (x) = J2 x n ^ + - +n f = 1 j- , (37) 

0<m<^<...<n f ^ - - 1 X 1 " X 3 ) • • • (1 - Xf) 

familiar from the theory of partitions (7j. In fact, we may employ the symmetry in the relabelled 
indices. To see how this works, let us symmetrize the case / = 2: 

0(ni < n 2 ) -> -{6{m < n 2 ) + 6(n 2 < n x )) = -(1 + 9{m = n 2 )) . (38) 

This can be obviously extended to larger /: the inequalities lead to a combination of terms with 
no restriction, terms where two labels are equated, terms where three labels are equated, terms 
where four labels are grouped in two pairs of equal ones, and so on. Using the function Z(x), 
we can convieniently determine the various coefficients by working out how Zf(x) can be split 
up in the corresponding way. Again for the case / = 2, this means writing 

Z 2 (x) = a 2 + - 2 , (39) 
(1 — x) z (1 — x z ) 

and solving this for general x gives, indeed, a = (3 = 1/2. As usual in the theory of partitions, a 
result for general / is prohibitively complicated, and therefore we give only the first few values 
of/: 

M X ) = 777 7777 + 7^ 



2(l-x) 2 2(l-x 2 ) ' 

_ , - 11 1 1 11 

Z ^ x ) = 7777 « + o7, 777 77 + 



3^ 



6(l-x) 3 2(l-x)(l-x 2 ) 3(l-x 
Z ^ x ) = 7777 777 + 777 77777 7777 + 7; 



24(l-x) 4 4(l-x) 2 (l-x 2 ) 3(l-x)(l-x 3 ) 



11 11 

+ -TZ oTo" + 



(1-x 2 ) 2 4 (1-x 4 ) ' 
11 1 1 1 1 

Z 5[ x ) = 77771 7777 + 7777 777777 77777 + 



120 (1-x) 5 12(l-x) 3 (l-x 2 ) 6(l-x) 2 (l-x 
11 11 11 



3^ 



4(l-x)(l-x 4 ) 8(l-x)(l-x 2 ) 2 6(l-x 2 )(l-x 3 ) 

■ (40) 
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The result for A{x) in these cases is therefore: 

A x {x) = e x H 2 (x) , 
A 2 (x) = y (h 2 {x) 2 + H 4 (xj) , 

A 3 (x) = j(H 2 (xf + 3H 2 (x)H A (x) + 2H 6 (x 

M{x) = ^(tf 2 (*) 4 + 6# 2 (x) 2 #4(x) + 8F 2 (x)F 6 (x) 
+3H A {xf + 6H 8 (x) 



A 5 (x) = —(h 2 (x) 5 + WH 2 (x) 3 H 4 (x) + 20H 2 (x) 2 H 6 (x)+30H 2 (x)H 8 (x) 

+15H 2 (x)H 4 (x) 2 + 2m i {x)H % {x) + 24H w (x)) . (41) 



These identities can easily be checked explicitly to modest order in x. Here, we have introduced 
the class of generalized hypergeometric functions 

/ T n\m 



n>0 



71! 



Obviously, H\{x) = e x , while H 2 (x) = Iq(2x). The rest of the cases are treated in a similar 
fashion. The number of dinstinct amplitudes that need to be calculated contains: 

B{x) = r (2x)A f _ x (x) (43) 
C(x) = A / _ 1 (x).{2/^(2x)-Jo(2x)} (44) 
V{x) = A f „ 2 (x)(l' (2x)) 2 (45) 

The generating function Q for the dinstict amplitudes can be written as in (|I3|): 

Q{x) = A f (x) + 2I (2x)Af-i(x) + 4B(x) + 2C{x) + &D(x) (46) 

Note the occurence of coefficients 4 and 6, which are dictated by a careful examination of the 
differing initial states that contribute. Some numbers for the case of / = 3, 4, 5 flavours, are 
shown in the following table 
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Total number of dinstinct amplitudes 


n 


/ = 3 


/ = 4 


/ = 5 


2 


35 


35 


35 


3 


123 


123 


123 


4 


777 


777 


777 


5 


3,853 


3,853 


3,853 


6 


25,327 


31,087 


31,087 


7 


139, 975 


200,455 


200,455 


8 


870, 485 


1,676,885 


1,999,445 



Note that for small n the numbers coincide: this is due to the fact that, n = 4, say, allows no 
room for 4 different quark flavours to occur in one diagram, and the difference between / = 3 
and / = 4 can therefore only appear for n > 6. This is reflected in the fact that Af(x) coincides 
with Af-i(x) up to the x 2 $ term. 

In order to estimate the large x expansion of the generating function, as in section 2, we 
note that Hf{x) is 



Hf(x) ~ 



a 2xf 



l + 0(~ 

X 



x — > oo 



(47) 



(4vr) 1 /2 x i/2 

In the Appendix we show how the asymptotic expansion can be computed systematically and 
using this, we can approximate Af(x) by 



A f (x) = j { Hl{x) + \ Ht 2 {x)Ht{x) 



(48) 



the previous equation, gives corrections to the leading result. Using this and the derivatives 



taking into account equations ( pl[) and keeping only the largest powers of x. The second term in 
the previous equation 
of the Bessel function 

l' Q (2x) , i#(2x)~ff 2 (z)(l + 0(i)) 



x — > oo 



(49) 



we can estimate the large x expansion of the functions in (43-45): 



B(x) 
C(x) 



A f -i(x)H 2 (x) = jj^-jyH£(x) + 2(/ l 3) , 
A f -i(x) (^'(2x) - I (2x)) = A f -i(x)H 3 (x) = B(x) 



H^ 2 (x)H 4 (x) 



V{x) ~ ^, 2 (x)(/;) 2 = 7T -L_g/( 3;) + 1 Ht 2 (x)H^x) 



(7-2)! 



2(/-4)!" 



and the generating function: 



1 + 2/ + 6/ 2 y 6/ 2 -22/ + 2 ir/ „ 2 ^ 



/! 



2(/-2)! 



(50) 

(51) 
(52) 

(53) 
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We can also compute the coefficients of the asymptotic expansion of Q{x). To this end we 
calculate the coefficient in the expansion of the functions H2(x), H^{x) using the Borel transform. 
In particular for H2(x) we define 

P(x) = e x H f 2 (x) = Y, K nx n (54) 

n 

To estimate the coefficients K n we perform a transform on P(x): 

/ dye-yyf/ 2 P(xy) = £ K n T(n + L + l)x n (55) 
Jo ^2 

and we get 

i (i + 2/r+// 2 

K n ~ /, Nf/o — ; 1 r (5o) 

(4tt)// 2 r(n+| + l) 

Similarily, for H^{x), approximated by 

H ^ X ) ~ (32^3)1/2^/2 > * ~> °° ( 5? ) 

we use Q(x) = e^-ff^ {x)H^{x) = ^ n L n x n . Performing a Borel transform we get 

/ dye-yy2 + 2Q(xy) = Y,L n Y{n+ J - + -)x n (58) 
Jo 2 2 

and for the coefficients 

1 (1 + 2/) n + //2 + l/2 ^1 + 2/V/ 2 T(n+^ + l) 

L ™ ~ /, w/o f 57 = K ™ o- 77 — — ? 57 ( 59 ) 



(47r)//2^fr(n + | + |) "I 2i j r(n + { + § 
Using these coefficients we can estimate the coefficients for the generating function in (53). 



Appendix: Asymptotic form of H m (x] 

Here we study the asymptotic form of the function H m (x), which was defined as 

mn 

= E M)^ ( 6 °) 

n>0 ^ 
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One can easily see that the following relation holds 

1 f dz x 
H m (x) = — j> —H rn -i(xz)Hi(-) 



1 f.. ( f^l...^zl Hl{xZl)Hl{xZ2) ... Hl{ E ) (6i) 

J J zi z m -i z 1 z 2 ---z m -i 



(2m) m ~ 1 

If we put Zi = the integral becomes 

1 t>2ir r27T 

H m {x) = j^tJ q #i-#™-ie^ (62) 

where 

W — g»01 _|_ . . . g«0m-l _|_ g-i(<^l+— <^m-l) (g3) 

We can estimate this integral by using the saddle point approximation. The first few derivatives 
of W are 



dW 

d<Pk 
d 2 W 
d(f)kd<f>t 
d 3 W 
d<j>kd4>ed(l)p 



l ^g^fc _ e ~ i (4>i+—4>m-i)^j 
- (f^^kt + e - iW,1+ -^ m - l) ) , 

{e i ^8 Mpq + e- i ^ + - ( l > ^) (64) 



d4>kdcj)ed(j) p d(f> q 

The saddle point can be found from the first derivative, and it is the solution of the equation 

e i<p + e -i(m-i)4> = _ e im<t, = 1 ^ (j)= *L k fe = 0, 1, . . . m — 1 (65) 

m 

The value of xW at the saddle point is x((m — \)e % ^ + e^) = mxe'f The saddle point that 
gives the largest real part of mxe 1 ^ is the one that dominates. We see that the function has an 
m-fold symmetry: if we restrict ourselves to \arg(x)\ < ^ the saddle point that dominates is 
= 0. The derivatives now take the values: 



d 2 W 

d(p k d(pe 
d 3 W 
d(p k d(f)ed(j) p 
d 4 W 
d(f)kd(j)ed(j)pd4> q 



- (Ski + 1) , 
-i (he P - 1) , 

(Skipq +1) , • • • (66) 
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and the exponent is 

xW = mx- {5kt + 1) 4>k<Pi ~ -7- Yl ( S M P ~ !) <Pk<p£(f> P + 7Ta ( S Mpq + !) <Pk4>£<Pp<Pq H 



(67) 

This is reminiscent of a zero-dimensional scalar field theory with vertices of arbitrary multiplicity, 
with the Feynman rules 



x V m 



^ (<W/3 + 1) 




[A a 

X 



ix (S^app - 1) 




where S^va = S^S^, S^aft = 8 ^5^5^, and so on. We can use the familiar tools of field 
theory to evaluate the integral. The first subleading term is computed by taking into account 
the two-loop diagrams that contribute. The result is 



5O— 0+|OO+n© = 



+ (m - l) 2 _ (m - l)(m - 2) = m 2 - 1 
8mx 12mx 24mx 



where the factors in front of the diagrams are symmetry factors. The next subleading term 
can be computed by including three-loop graphs. Due to the fact that — Q = 0, there are 
8 non-zero connected three-loop diagrams, and in addition to these we must also include the 
disconnected diagrams that are shown below. The result for the next term in the expansion of 
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H m (x) (including the symmetry factors shown below) is: 



+ 48 



+ s®4O0^(oo + e;" 



m-l)(m 3 + 289m 2 -1129m + 1175) (69) 



1152m 2 x 2 

The result for the asymptotic expansion of H m (x) to this order is: 

TfhX 2 1 1 1 

fl" m (x) ~ - 6 = {l+ m ~ + ^(m-l)(m 3 + 289m 2 -1129m+1175) + C(^)} 

V ^ Vm(2vrx) m - 1 24mx 1152m 2 x 2V A 7 x 3 " 

(70) 

and higher terms can be obtained in a similar way. 
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